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Abstract 

We show that the Lindstedt-Poincare perturbation theory is always a rehable technique in the 
region of small coupling constant. The harmonic balance result, on the other hand, if expanded in 
the perturbation parameter may lead to incorrect results. 
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I. INTRODUCTION 



A senior undergraduate course in nonlinear dynamicsii^ abounds in various kinds of ap- 
proximation procedures. There is the Lindstedt Poincare method, the harmonic balance, 
the multiple time scale analysis, Bogoluibov-Krylov technique and so on. It is interesting 
to note that for the cubic anharmonic oscillator the Lindstedt-Poincare perturbation theory 
answer given in the text by Landau and Lifshitz differs from the harmonic balance answer 
of Jordan and Smith if the Jordan and Smith result is expanded in a perturbation series in 
the cubic nonlinearity. In view of this difference we re-examine the Lindstedt-Poincare and 
the harmonic balance technique and compare with the perturbative expansion of the exact 
integral for the frequency of the motion. The Lindstedt-Poincare expansion is in agreement 
with the perturbative expansion of the integral. We point out how the harmonic balance 
method can be combined with perturbation theory to yield an answer in agreement with 
the Lindstedt-Poincare expansion. As another example of a possible discrepancy, we con- 
sider the periodic orbit of the Lotka-Volterra model. Application of the Lindstedt-Poincare 
technique yields an answer which differs from the perturbation expansion of the harmonic 
balance result given in Jordan and Smith. A numerical determination of the frequency shows 
agreement with the Lindstedt-Poincare method. It is apparent that the Lindstedt-Poincare 
method is always reliable. A perturbative expansion of a non-systematic harmonic balance 
answer may not give the right perturbation series when the equation of motion does not have 
reflection symmetry. To get the low coupling constant part right, harmonic balance has to 
be used in conjunction with perturbation theory. In case there is reflection symmetry, the 
frequency, determined by harmonic balance is automatically in agreement with perturbation 
theory. We point out that knowing the correct perturbation theory is vital to constructing 
a globally valid formula for the frequency in an equivalent linearisation technique. 

II. EXACT INTEGRAL AND LINDSTEDT-POINCARE PERTURBATION THE- 
ORY 

We use the example of the anharmonic oscillator 

X + uj'^x + Ax^ = (1) 
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to illustrate Lindstedt-Poincare, harmonic balance and the equivalent linearization tech- 
niques. This system is in principle exactly solvable as we can write the integral of motions 

as 

- I — I + -uj X + —X = constant [2) 

The constant can be evaluated by noting that i; = at the turning point i.e. at x = a, 
where a is the amplitude of motion. This leads to 

(|)'=c.^(a^-x=)+^(a'-x') (3) 

and consequently 

dt = — — (4) 

[uJ^a^ - x^) + ^ {a'' - x%^ 

For the symmetric potential, that we have here V{x) = ^cu'^x'^ + jx'^, the time period T is 
found from Eq. (jlj) as 

T r dx 



1/2 

de 



4 Jo u{a^-x^y^^ [1 + 2^(a2 + x2 



(5) 



/o ^[l + ^(l + sin2^)] 

We can express the right hand side of the above equation in terms of the elliptic function 
For jluP' <^ 1, we can expand the integrand as, 



T r''^ de /-"/^ de 

3Aa2 



/ / — l + sin^e) + ... 

Jo ^ Jo ^ 

(6) 



TT 

2^ 



1 - 



while for Xa /lu 3> 1 



T 1 r^^ de 



4 uJn V Aa2 (l + 3in2^^)l/2 



^ 7r(i/4) (7) 



The frequency VL of the oscillations is = 27r/T and is found from Eqs. ([6]) and ([7]) to be 
given by 

n = uj + - — + ... (8 

8 uj^ 
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and 



r(i) 



UJ 



The Lindstedt-Poincare perturbation theory for Eq. (JT]) proceeds by expanding 

X = Xq + Xxi + A^X2 + . . . (10) 
and writing the actual frequency Q as 

= uj^ + Xul + x^ujI + . . . (11) 

In terms of the above expansions, Eq. ([T]) can be written as 

io + Ail + A^i2 H h (xq + Axi + A^X2 + . . . ) 

= -A (xo + Axi + . . . ) + (xo + Axi + . . . ) + A^cj^ (xq + Axi + . . . ) (12) 

Equating the same powers of A on either sides, we get 

A°: Xo + n^Xo = (13a) 

A^: xi + fi^xi = + w^xo (13b) 
A^: X2 + f2^X2 = — 3xqXi + tu^Xi + u;|xo (13c) 
The solution of Eq. fll3al) for the initial conditions Xq = ^Jo = at t = is 

xo = Aq cos fit (14) 

With Xo obtained, Eq. (]13bll becomes 

Xi + fi^xi = — cos^ fit + ciJ^Ao cos fit 

= [ ^ + cu?Ao j COS fit - ^ COS 3fit (15) 

The drive with frequency fi on the right hand side causes a spurious resonance in the system. 
This would cause x\ to diverge. In order to have a finite xi, we need to remove the resonance 
causing term from the right hand side of Eq. fll5l) . This is done by the choice 

ul = ^ (16) 

The equation for xi can now be solved and for the frequency fi we have 

^2 ^ ^2 ^ ^ _ _ _ (17) 



or 

n = u + '-^^ + ... (18) 

o OJ 

which is the same as that shown in Eq. (jH]) 

Harmonic balance, on the other hand, requires the expansion (without loss of generality, 
we consider a solution with x = 0) at t = 

X = ao + ai cos Qt + ^2 cos 2Qt + ^3 cos 3Qt + . . . (19) 

Inserting this solution in Eq. ([1]) 

— ri^ai COS ^It — 4fi^a2 cos 2i^t — Qfi^as cos Sfit + . . . 
+00"^ (ao + ai cos Qt + 0:2 cos 2Qt + . . . ) 

+A (ao + 3aiao cos fit + Sa^ao cos^ fit + af cos^ fit + ... ) =0 (20) 

Equating the coefficient of each harmonic separately to zero: 

ao = (21) 
]a 



= uj^ + Ma (22) 



III. THE ANHARMONIC OSCILLATOR 

We begin this section with the cubic oscillator 

x + uj^x + Xx^ = (23) 

The potential is V{x) = \uj'^x'^ + ^x^ and we note that while it is unbounded for negative 
x(A > 0), it has positive maximum at a; = — cu^/A and if the total energy E (determined by 
initial conditions) is such that it is less than A, then there will be periodic motion with 
frequency fi, with the time period given by 

r dx 

T = / , = (24) 

where ai and 02 are the consecutive negative and positive real roots of the amplitude equa- 
tions 

cj^a^ Aa^ 

E ■ 
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or a' = — - -— 25 



For a perturbative expansion of the integral, we first need to find ai and 02 from a pertur- 
bative determination of the roots of the cubic of Eq. fl25|) . Expanding the root A as 



A = A0 + XA1 + AM2 + . . . 



(26) 



we have from Eq. (125!) 

Al + 2XA1A2 + X^Al + 2X^AoA2 + 
Equating equal powers of A on either sides, 



^0 



(27) 



Al 
A2 



2E 

^0 
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(28a) 

(28b) 
(28c) 



At the zeroeth order ai = —\J^ and 02 = ^J^■ To (9(A^), we have from above 



ai 



a2 



2E ^2E 3^2 5 1 f2Ey^^ ^ 



^3cj4 ^ 18a;H^V 



2E ,2E ,0 5 1 /2E\ 



3/2 



- A— + A 



3cu4 18 uj^\uj'^ J 
Noting that ai < and 02 > 0, we can split the integral in Eq. fl2^ as 



(29) 
(30) 



T 
2" 



+ 



ail 



'ail 



+ 



02 



2E - uj'^x^ + ^ -^0 J2E - cu2x2 - 2^ 



dx 



dx 



^cu2 (a2-x2)-f (|ai|3-x3) -^0 

1 r^/^ 1 r''^ 



2 A|ai / 1 



3 oj^ V cos 



u;2 (a2 — x2 



^ (a^, - x^ 



3 V'^2 



I _j_ 2 Aa2 /' l~sin^ ' 



1 /"'r/2 

a; 



3 o;-^ V cos^ 
2 



sm'^ 6^ ^ 



cos2 6 



+ 1- 



Aa^/l-sin^^A 3 / 2A A V 1 - sin^ ^ 



3u;2 \ cos2 9 



+ 7: 



8 V3cj2 



cos2 9 



(31) 
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Noting that |ai| and 02 differ at 0(A), we see immediately that the corrections to the leading 
order answer of tt for the integral comes at O(A^) and is given by 



2~ uj uj 9uj^ Jo 



1 -sin^e 3 /I -sin^e' ^ 



cos^ 6 ~'~ 2 V cos^ 6 



de 



TT A^ AE flh \ TT X^5 E 

- H a TT = - H rTT 

UJ UJ 9uj^ \ 8 J UJ UJ Quj^ 



(32) 



To implement the Lindstedt-Poincare scheme, we expand as before(i7 is the real frequency 
of oscillations) 



X = Xq + Xxi + \^X2 + . . . 

= uj'^ + XujI + X^ujI + ... 
Inserting in Eq. (l23ll and equating identical powers of A on either sides 



A°: 
A^: 
A^: 



Xq + VL^Xq = 

xi + ri^xi = 

X2 + fl'^X2 = —2XoXi + Ujfxi + UJ^Xo 



2 I 2 

-Xq + UJ-^^Xq 



The solution of Eq. (l36i) (initial conditions a;o(0) = 0, xq{0) = 0) is. 



Using this in the Eq. (j37ll we have. 



xo = A cos Qt 



^2 



xi + fi^xi = — — (1 + COS 2nt) + ujf cos rit 



(33) 
(34) 



(35) 
(36) 
(37) 



(38) 



(39) 



As before, there is a resonating term on the right hand side and removal of this requires 



t^i = 



(40) 



The solution xi is now 



Xi = Bi cos Vtt + B2 sin Vtt 



A^ 



cos 2VLt 



2n^ 6fi2 

Using xi = xi = at t = 0, we get B2 = and Bi = With this information, Eq. ([3 
reads 



(41) 



A^ 

X2 + f2^X2 = 777^ (1 + COS 2Qt) - — - (cos 3Qt + cos Qt) + uj\A cos (42) 



31^2 



Once again the resonating terms need to be removed from the right hand side and this leads 

to 

-5 - («) 

The frequency thus to O(A^) is = lu — + • • • agreement with Eq. (15^ and the 

text of Landau and Lifshitz. 

A harmonic balance approach requires an expansion 

X = ai + a2 cos Qt + cos 2Qt + . . . (44) 

The equation of motion becomes 

— fi^a2 cos ^It — 4i7^a3 cos 2^lt + ■ ■ ■ + u;^ (ai + a2 cos 2^lt) 

+A (^al + 2aia2 cos + y (1 + cos 2fit) + 2aia3 cos 2^^ + 02^3 cos (^t + . . . ^ =0 (45) 

Setting the coefficient of each harmonic separately to zero, we get 

2 9 
Y + «? = -y«i (46) 

a2 {uj^ -n^ + 2Xai + Xas) = (47) 
as (c^' - 4fi') + ^ + 2Aaia3 = (48) 

If we ignore as, as recommended in Jordan and Smith, then the perturbative solutions yields 

- («) 

\ 2 2 

and ^ + 2Aai = - ^ (50) 

in agreement with the result quoted there and in disagreement with the (9(A^) term of 
Landau and Lifshitz and the direct result obtained in Eq. (|32l) . Keeping the as term, 

Aaa Aai ^/,2\ 

Using this in Eq. (I48l) we finally get 

l^^ = c^^-^^ (52) 
6 to'"' 

which is the correct answer at (9(A^). 
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IV. LOTKA-VOLTERRA MODEL 



The predator-prey system is a two variable dynamical system with x standing for the 
number of prey at time t and y the number of predators. Left to themselves, the number 
of prey grows (rabbit feeding on grass) while the number of predators decreases (foxes left 
with no food). An interaction between the two species causes the number of prey to fall and 
predators to increase. An interaction between the two species causes the number of preys 
to fall and the predators to increase. The interaction effect is jointly proportional to the 
number of prey and predators. These considerations led to the Lotka-Volterra model 

dx 

— = X — xy 

t (») 

The model has two fixed points x = y = Q which is unstable and x = y = Iwhich is a 
center, i.e. small deviations from it execute a periodic motion. In reality there is a periodic 
orbit about the point (1, 1) and it is the time period of that orbit which we would like to 
determine. 

The first step is to introduce variable Xi, X2 centered around (1, 1) where Xi = x — 1 and 
X2 = y — 1 and in terms of these variables, the equation of motion become 

Xi = -X2 - X1X2 

(54) 

X2 = Xi+ X1X2 

It is clear from the above that for small Xi,X2 the linearized system x'l = —X2 and X2 = Xi 
corresponds to simple harmonic motion with frequency unity. The nonlinear terms make the 
time period, amplitude dependent. This is what we set out to find. To do so we introduce 
a parameter A in front of the nonlinear term and the write Eq. fl5^ as 



Xi = —X2 — \X1X2 

(55) 

^2 = Xi + XxiX2 

For A = the dynamics is that of a simple harmonic oscillator with xi = Aicost, X2 = 
A2Cost where Ai and A2 are the initial values of xi and X2. For A 7^ 0, we explore the 
possibility of existence of this periodic solution. The frequency of motion will no longer be 
unity but some frequency Q which needs to be determine. We anticipate for small A 

n = l + XuJi + \^UJ2 + ■■■ (56) 



while 

2^1,2 = a;io,2o + Aa;ii,2i + A^xi2,22 + • • • (57) 

We rewrite the Eq. fl55l) as 

Xi = —QX2 — \X1X2 + (XuJi + + . . .) X2 

(58) 

X2 = fiXi + AxiX2 + (AcJi — A^Co'2 + ■ . .) X2 

We now expand Xi 2 as in Eq. (157|) and insert the series in Eq. (l58l) . Collecting and equating 
identical powers of A on either side from the resulting equation, we get 



A° : Xio + ^X2o = 

X20 - ^Xio = 

A^ : ill + ^X2i = -X10X20 + 1^13^20 

X21 — ^Xii = 2:10X20 — ^^iXlQ 

A^ : X12 + ^X22 = -2:10X21 - a:2o2;ii + cl'iX2i + 1^22:20 

X22 ~ ^Xi2 = a;ioX2i + X2o2;ii — iOiXii — ct;2a;io 



(59) 



(60) 



(61) 



Clearly at (9(A°) the solution is, xio = AicosQt + y42 sinfit, X20 = Ai sin fit — y42Cosfit 
where Ai and A2 are taken to be the initial values of xi and X2- At 0{X^), we have 

_ A1A2 

Xii + fix2i = — - sin 2ilt H — cos 2ilt + uji {Ai sin fit — A2 cos fit) (62a) 

_ j^2 A^Ao 

X21 - fix2i = — - sin 2fit cos 2fit + ooi {Ai cos fit + A2 sin fit) (62b) 

The homogeneous equation has the solution of the form A cos fit + i? sin fit which shows 
that the cos fit and sin fit in Eq. fp2]) will cause the system to resonate, removal of this term 
requires uJi = 0. Now from Eq. fl62l) 

xii + fi^xii = fi {A1A2 - {Al - Al)) cos2fit - fi (2A1A2 + ^{Al- Al)^ sin2fit (63) 

yielding 

xii = A cos fit+S sin fit — ^ [A1A2 - {Al - Al)] cos2fit+^ [{Aj - Al) + 4A1A2] sin2fit 

(64) 
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At t = 0, xii = which sets A = ^ [^1^2 - (Aj - Aj)]. From Eq. ([62]) we have, 

Qx2i = -BQcosQt ^ [A1A2 - {AI - Al)] sin fit 

-— [AA1A2 + (A^ - A^)] COS Qt - — - sin fit + cos fit (65) 

At t = 0, xii = leads to 5 = -3^^ [A1A2 + {Aj - Aj)]. So Finally we have 

xn = A cos fit + 5 sin fit + C cos 2fit + D sin 2fit (66a) 

X21 = A cos fit - S sin fit + cos 2fit + F sin 2fit (66b) 

where, 

A= ^[A,A2-{Al-Al)], B= -^[A,A2+{Al-Al)] 
C= ~^[A,A2~{Al-Al)], D= ^[^A,A2+{Al-Al)] (67) 
F= ^[-^A,A2+{Al-Al)], E= -^[AA2+{Ai-Al)] 
Now at (9(A^) we have the system of equations 



xi2 + fia;22= -(a;ioa;2i + a;2oa;ii) +u;2a;2o 
i;22 - ^xi2 = (a;ioX2i + X2oXu) - c^2a;io 



(68) 



We first compute the term (xioa;2i + a;2oa;ii)- Using results obtained for xio, X20, ,a;iiX2i we 
get 



3^10^21 + X20XU = {Al COS fit + A2 sin fit) {A cos fit + -B sin fit + C cos 2fit + D sin 2fit) 
+ {Al sin fit - A2 cos fit) {A cos VLt - B sin fit + E cos 2fit + F sin 2fit) 
= [higher harmonics] + ^ {AiE - A2C + AiD + A2F) cos fit 

+ ^ {A2E + AiC - AiF + ^2^) sin fit (69) 

We need not pay heed to coefficients of the higher harmonics as we just need to look at 
terms that will cause the system to resonate. Such terms must be put to zero from physical 
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considerations. From Eq. (l68l) we get 



xi2 + = {xiqX2i + a;2oa;ii) + Vtu2XiQ - (xioX2i + a;2oa:ii) + u;2a:2o 



d 

(X10X21 + a;2oa;iij + ilu2XiQ - 
= {AiE - A2C + AiD + A2F + A2E + AiC - AiF + A2D) cos l^t 
(^2^ + AiC - AiF + A2D - AiE + A2C - AiD - A2F) sinl^t 
+2Vtuj2 {Ai cos cut + A2 sin l^t) + ^ {AiE - AaC + + A2F) sin fit (70) 

Now to get rid of the resonance creating terms cos VLt and sin VLt we must take their coeffi- 
cients to be zero. Thus we get 

2UJ2A1 = i [-AiE + A2C -AiD- A2F - A2E ~AiC + AiF- A2D] (71) 

2(^2^2 = I [-A2E -AiC + AiF- A2D + A^E - A2C + A^D + A2F] (72) 
From Eqs. (ITTjl and (172!) we obtain 

2uj2 (Ai + A2) = (F - C) - A2 [E + D) (73) 

Now using the expressions for Ai, ^2, C, i?, and F we get 

2072 (Ai + A2) = (Ai - A2) - ^ (Ai + A2) (74) 



Finally we have the expression for 102 



^1 + ^2 



The perturbative result for fi up to O(A^), after setting A = 1 is 



\A 



2 



where A is the amplitude of the limit cycle. If xq and be the initial values of x and 
then we can write 

^ ^ 1 _ (^o^l)!_H(yo^ (77) 

We have checked this result numerically. The results are shown in Fig[l]. The good 
agreement between the computed frequencies and obtained from Eq. ( |76l) is apparent. 
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FIG. 1: Initial number of population density modifies its periodicity. Initial potulation density 
(x=y) is plotted with corresponding frequency (uj) of its oscillation. Dotted line: Theoretical 
prediction. Solid Line: Numerical result 

V. EQUIVALENT LINEARIZATION AND CONCLUSION 

In this section, we conclude by explaining the equivalent linearization process which spans 
the entire range from low to high expansion parameter. As will be apparent in the immensely 
useful process the input coming from the low order perturbation theory is vital. 

The technique of equivalent linearization proceeds by replacing the non-linear term 
by a linear term a{x'^)x, where a is a number of 0{1) and (x^) is the average of x'^ over 
one cycle. We now simply write Eq. ([1]) as the equation of motion of a simple harmonic 
oscillator 

i + n^x = where fi^ ^ uj^\a{x^) (78) 

The first question is how do we fix the value of a? In terms of the potential, we have replaced 
jx'^ by and if we demand that this statement be true at least on an average, then 

we would require 

(x^) = 2a{xY (79) 

and using periodic solution x = Aq cos nt, this leads to a = 3/4. Writing (x^) = ^, Eq. ([78]) 
shows = c<j^ + IAAq or 

n = u;(l + 4a41 (80) 



16 Cj2^ 

which differ from the Lindstedt-Poincare answer at 0{X) by a factor of two. In using the 
equivalent linearization technique, it is best to use it in conjunction with the Poincare 
Lindstedt perturbation theory to fix the parameters of the linearization scheme. In this 
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particular case the agreement with the correct perturbation theory requires fixing of a as a = 
3/2. The strength of the equivalent linearization lies elsewhere. With the equivalent simple 
harmonic oscillator of Eq. (175]) the integral of motion is + ^Vt^x^ which is recognized as 
the energy of the oscillator. The energy is a physical quantity which is fixed by the initial 
condition and the primary quantity for the dynamics as opposed to the amplitude. For 
the simple harmonic oscillator ^fi^A^ = where A is the amplitude while for the original 
oscillator E = \uj'^a^ + with a being the corresponding amplitude. Returning to our 
equivalent oscillator now, we write {x^) = ^ = ^ and using Eq. (ITHl) with a = 3/2, find 

"^--^ + 1^ (81) 



with the result 

2 



UJ 



2 + Vuj^ + QXE 



(82) 



Expanding in powers of A, this gives the correct answer to (9(A) as it should but goes further 
than any perturbation theory by giving a meaningful large A limit. For XE/uj^ ^ 1, Eq. (!82ll 
shows 



1/4 



n^^^j iXEf' (83) 

which is to be compared with the exact answer of ~ [V27r/r(i)] (AE)^/^ The two 
prefactors differ by a meager 10%. The entire range of A can be handled by Eq. ( l82l) ). 
wherein lies the strength of equivalent linearization. It is important to note that Lindstedt- 
Poincare technique provides a necessary input for the success of equivalent linearization. 
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